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Abstract. The most puzzling aspect of the radio emission from LS I +61° 303 is that the 
large periodic radio outburst, with period equal to the orbital one, occurs very displaced 
from periastron passage, nearly at apoastron. In 1992, Taylor, one of the discoverers of this 
source, together with his collaborators proposed a model of a compact object in an eccen- 
tric orbit accreting from the equa torial wind of the Be star primary. The application of this 
model bv lMartf & Paredes \ 19951) predicts one ejection at periastron and a second more dis- 
placed ejection along the orbit. The first ejection should correspond to weak radio emission, 
because of strong inverse Compton losses of the emitting electrons due to the proximity 
to the hot Be star, whereas the second ejection, quite displaced from the star, would corre- 
spond to a strong radio outburst, that one indeed observed. Corroborated along the years by 
numerical computations, simulations and gamma-ray observations, until now this two-peak 
model could not be proved in the radio band, because of the negligible emission around 
periastron. We show here, that the radio spectral index based on the ratio of flux densities 
is the unique tool to monitor activity of LS I +61°303 in the radio band around periastron. 
The analysis of the radio spectral index over almost 7 years of Green Bank Interferometer 
data results in a clear double-peaked spectral index curve along the orbit. This result gives 
finally observational support at radiowavelengths to the two-peak accretion/ejection model 
for LS I +61°303. Moreover, the here shown comparison of the two-peak curves - the ra- 
dio spectral index curve and the Fermi-LAT gamma-ray curve - indicates a new interesting 
hypothesis on the electron population responsible for the gamma-ray emission. 

Key words. Radio continuum: stars - X-rays: binaries - gamma-rays: observations - X- 
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1. Introduction 

The TeV-emitting source LS I +61°303 is a 
X-ray binary system where a compact ob- 
ject travels through the dense equatorial wind 
of a Be star. The most typical peculiarity 
of LS I +61°303 is a large periodic ra- 
dio outburst toward apoastron. In Fig. 1 we 
see 6.7yr of Green Bank Interferometer data 
folded with the orbital period P\ =26.496 d 



(<£: with to=JD 2443366.775) dGregorvl 

2002). Periastron passa ge corresponds to or- 
bital phase <5=0.23 fo r ICasares et all (120051) 
and to <D = 0.275 for lAragona et al. (2009). 
As one can see in Fig. 1, the outburst is clus- 
tered around i> ^0.6, i.e. almost apoastron. 
The broad shape of the light curve of Fig. 1 is 
due to variations of the orbital phase and am- 
plitude of the outburst, both changing with a 
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Fig.l. Radio light curve of LS I +61303 (6.7 yr 
of Green Bank Interferometer data ) vs orbital phase 
(M as si & Kaufman Bernadd 2009) 



long-term period P 2 =1667 d (lGregorvll2002h 
(0: £=^). IZamanov & Martil (120001) demon- 
strated a modulation on the same timescale in 
the H a emission line. This latter result strongly 
suggests that the long-term modulation in the 
radio emission is related to changes in the Be 
star equatorial disk pro perties. In particular, 
iGregorv & Neishl d2002l) suggest that the long- 
term modulation in radio properties may stem 
from periodic ejections of a shell (density en- 
hancement) of gas in the equatorial disk of the 
Be star. 



2. Two Peak Accretion Model 

One of the fundamental questions concerning 
the strong periodic radio outbursts at O ^0.6 
of LS I +61°303 has been: why are they shifted 
with respect to the periastron passage (i.e. O = 
0.230-0.275)? 

The accretion rate M oc £^£± > where p w i n( j 

l 'rd 

is the density of the Be star wind and v re i is 
the relative speed between the accretor and the 
wind, is in fact propo rtional to the dens ity of 
the accreted material dTavlor et al.lll992h . The 
highest density is obviously at the periastron. 
The explanation for the shift is that the orbit 
of LS I +61°303 is quite eccentric, e=0.54 - 
0.7 dAragona et al.ll2009t ICasares et al.ll2005l) . 
In an eccentric orbit the different relation- 
ship for density and velocity (inversely pro- 
portional to the power of 3) creates two peaks 
in the accretion rate curve, one at periastron 



because of the highest density, and a sec- 
ond one when the drop in density is compen- 
sated by the decrease in v elocity towards apas- 
tron. [javlor et al] (1 19921) computed the accre- 
tion rate curve for different eccentricities and 
showed that two peaks begin to appear for an 
eccentricity above 0.4. Whereas the first peak 
is always toward periastron, the orbital occur- 
rence of the second accretion peak depends 
on variations of the wind of the Be star with 
the period P 2 =1667 d. lMartf & Paredesld 19951) 
computed the accretion rate curve for different 
wind velocities and showed that for a stellar 
wind velocity of 20 km/sec the two peaks be- 
come rather close to each other, whereas for 
a wind velocity of 5 km/sec they have an or- 
bital offset of AO = 0.4. For periastron at 
O 0.2-0.3 the second peak may therefore 
appear at O ~ 0.6 - 0.7, as indeed occurs. 
iMartf & Paredesl (1 19951) have shown that both 
peaks are above the Eddington limit and there- 
fore one expects that m atter is ejected twice 
within the 26.5 d interval. Ikomero et al.l d2007l) 
applied a smoothed particle hydrodynamics 
code to develop three-dimensional, dynamical 
simulations for LS I +61°303 and found that 
indeed the accretion rate has two peaks per or- 
bit, i.e. a narrow peak at periastron, and a broad 
peak that lags the periastron passage by about 
0.3 in phase. 

Models and simulations predict therefore 
two peaks. Wh y do we observe onl y the 
second ejection? Marti & Paredesl d 1995b pre- 
dicted that near periastron the ejected rel- 
ativistic electrons are embedded in such a 
strong UV-radiation field that they loose their 
energy by the inverse Compton (EIC) pro- 
cess: no radio emiss ion but high energ y emis- 
sion is predicted. iBosch-Ramon et all d2006l) 
computed the inverse compton losses and 
the related light curves of emission for 
LS I +61°303 in the radio band a nd at high en- 
ergy. Bosch-Ramon et al ] d2006l) fixed the Be 
wind parameter to have the second accretion 
peak at O=0.5. Their results, are here shown in 
Fig. 2. In perfect correspondence to the accre- 
tion rate curve the radio light curve shows two 
peaks, one large outburst at O=0.5 and another 
smaller outburst at periastron. The high energy 
light curve shows exactly the contrary situa- 
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Fig. 2. Computed accretion r ate and light curves 
by iBosch-Ramon et alj d2006l> (from their Fig. 2). 
Bosch-Ramon et al. fixed the Be wind parameter to 
have the second accretion peak at 0=0.5. 



tion, again two peaks one at periastron and the 
other at 0=0.5, but the dominant peak is at pe- 
riastron. During the second accretion peak the 
compact object is at larger distance from the 
Be star and therefore inverse Compton losses 
are lower: the associated gamma-ray outburst 
is weaker and the electrons are able to emit 
stronger synchrotron radiation producing the 
larger radio outburst. 

Gamma-ray obse rvations confirm 
Marti & Paredesl dl995|) p redictions and 
(2006) calculations, 
detected by EGRET 
As discussed in 
and 
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Bosch-Ramon 
LS I +61°3 3 was 
dTavani et all Il998l) 
Massi (2004), Massi et al] d2005l) 



Massi & Kaufman Bernadd d2009h these data 
well support the hypothesis of a high-energy 
outburst at periastron: EGRET observations at 
0=0.18 during a well sampled full orbit show 
a clear peak at periastron passage. EGRET 
observations at 0=0.41, along with an increase 
of the emission again near periastron show 
ev en a second peak a t O ^ 0.5 (see Fig. 3 
in iMassietai] d2005[)'). Fermi-LAT observa- 
tions dAbdo et al.l l2009l) were performed at 
= 0.788 - 0.927. The Fermi light curve is 
characterized by a broad peak after periastron 
as well as a smaller peak just before apastron. 
Therefore, the first gamma-ray peak seems 
indeed to be persistent, whereas the second 
gamma-ray peak seems as predicted, to change 
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Fig. 3. LS I +6T303. Spectral index and flux den- 
sity of Green Bank Interferometer data at 8.3 GHz 
and 2.2 GHz vs orbital phase, <J>, in the interval 
0=0.0-0. 1. Note that the evolution from an optically 
thick to an optically thin spectrum occurs twice, 

iving the a vs <t> curve a doub le-peaked shape. 

Massi & Kaufman Bernadd 2009) 



as function of the long term variations of the 
Be star. 



3. Radio spectral index analysis 

From the previous section we see as the 
two-peak accretion model, including energetic 
losses, predicts at periastron a high energy out- 
burst due to IC along with a small radio out- 
burst, and associated to the second accretion 
peak a large radio outburst and possibly, de- 
pending on the stellar distance, a gamma-ray 
peak. The large radio outburst should follow 
the typical characteristics of microquasars: op- 
tically thick emission, i.e. a > (with flux 
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density S oc y"), follo wed by an opticall y 
thin outburst, i.e. a < dFender et al.ll2004l) . 
In microquasars the first type of emission, the 
optically thick radio emission, is related to a 
steady, low velocity, conical jet centered on 
the system. The following optically thin out- 
burst, called "transient jet", is due to shocks 
caused by the travelling of new highly relativis- 
tic plasma, generated by a transient, through 
the underlying, slow, steady flow (Fender et 
al. 2004; review by Massi in this volume). In 
Fig. 3 we see the spectral index and flux den- 
sity, at 8.3 GHz and 2.2 GHz vs orbital phase, 
0>, f or the GBI data in the interval 0= 0.0-0. 1 
(see lMassi & Kaufman Bernaddld2009l) for de- 
tails). Clearly associated with the large out- 
burst of LS I +61°303 is the predicted evolu- 
tion for microquasars, from an optically thick 
to an optically thin spectrum. At the bottom of 
Fig. 3, the light curve at 2.2 GHz reveals that 
the large outbust at cD ~ 0.8 is preceeded by 
another outburst at O ~ 0.7. The spectral in- 
dex at the top of Fig. 3, shows the different na- 
ture of the two outbusts: the minor outburst at 
<t> = 0.7 is optically thick, whereas the larger 
peak around <t> = 0.8 is an optically thin out- 
burst. In particular, the outburst at <S> = 0.7 
is related to an optically-thick-emission inter- 
val creating a broad peak in the spectral index 
curve. In the context of microquasars this inter- 
val corresponds to the emission from a steady, 
low velocity conical jet. Around periastron in 
the light curves at 2.2 and 8.3 GHz one sees 
only very small, barely detectable outbursts at 
<t> ~ 0.3 - 0.4. On the basis of the light curves 
alone one would never be able to associate 
these negligible peaks to th e predicted small 
radio o utbu rst at periastron bvlMartf & Paredei 
(1 19951) and lBosch-Ramon et all (120061) . When 
one, however, analyses the spectral index curve 
at the top of Fig. 3, one sees that O ~ 0.3 - 0.4 
indeed corresponds again to a broad peak in 
the spectral index curve. The evolution from 
an optically thick to an optically thin spectrum 
occurs clearly twice giving the a vs d) curve a 
double-peaked shape, as expected from a two- 
peak accretion curve. 

As shown in Fig; 5 of 

iMassi & Kaufman Bernadol (120091) this 
shape is not constant but changes during 



the long 1667 d cycle in agreement with 
iMartf & Paredesl {1995) computations for 
variable parameters of the wind of the Be 
star. This variation of the accretion rate 
curve/spectral index curve with implies that 
one can compare data of different epochs only 
when observed in the same interval. As an 
example one sees that at the top of Fig. 3, 
for = 0.0 - 0.1, the spectral index curve 
at O = 0.5 gives a < (i.e. corresponding 
to a transient jet); at the top of Fig. 4 for 
= 0.788 - 0.927, one sees that at the 
same orbital phase <t> = 0.5, it results a > 
(i.e. corresponding to a steady slow conical 
outflow). We compare here in Fig. 4 the 
Fermi-LAT gamma-ray curve with the spectral 
index for GBI data at other epochs than 
Fermi-LAT observations but in the same phase 
= 0.788 - 0.927. It is worth noting that both 
Fermi-LAT gamma-ray peaks correspond to 
intervals where the radio emission is optically 
thin. This correspondence would imply that the 
electrons responsible for the inverse Compton 
process creating the gamma-ray emission are 
those of the fast transient jet and not those of 
the slow outflow. 

4. Conclusions 

We analyzed the radio spectral index using 
6.7 years of GBI radio data of LS I +61°303. 
Our main conclusions are that the peri- 
odic (Pi = 26. 5d) large radio outburst of 
LS I +61°303 consists of two successive out- 
bursts, one optically thick and the other opti- 
cally thin. This creates one peak in the spec- 
tral index curve. In microquasars, the optically 
thick emission is associated to a steady jet and 
the optically thin outburst to the transient jet. 
We observe that along the 26.5 d orbit the evo- 
lution from an optically thick to an optically 
thin spectrum occurs also around periastron 
creating a second peak in the spectral index 
curve, giving the avsO curve a double-peaked 
shape. This result agrees with the predic- 
tions of the two peak accretion/ejection model, 
with the results of three-dimensional dynam- 
ical simulations and finally with gamma-ray 
data. All these results indicate a scenario with 
a first ejection around the periastron passage 
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Fig. 4. Radio and gamma-ray data observed at dif- 
ferent epochs but in the same interval ©=0.788- 
0.927 of the 1667 d periodicity. Top: Spectral index 
of Green Bank Interferometer data at 8 .3 GHz and 
2.2 GH z. Bottom: Fermi-LAT data from lAbdo et alJ 
(2009) 



with low radio emission, but high energy emis- 
sion, due to inverse Compton losses caused by 
the proximity of the B0 star, and a second ejec- 
tion far away from the Be star, with negligible 
losses and therefore with a well observable ra- 
dio outburst. Finally, the here given compar- 
ison of radio (GBI) and gamma-ray (Fermi- 
LAT) data seems to imply that the electrons re- 
sponsible for the inverse Compton process cre- 
ating the gamma-ray emission are those of the 
fast transient jet and not those of the slow out- 
flow. 
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